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ABSTRACT 

A number of recent investigations have revealed that transverse waves are ubiquitous in the solar chromo- 
sphere. The vast majority of these have been reported in limb spicules and active region fibrils. We investigate 
long-lived, quiet Sun, on-disk features such as chromospheric mottles (jet-like features located at the bound- 
aries of supergranular cells) and their transverse motions. The observations were obtained with the Rapid 
Oscillations in the Solar Atmosphere (ROSA) instrument at the Dunn Solar Telescope. The dataset comprises 
simultaneous imaging in the Ha core, Ca II K, and G band of an on-disk quiet Sun region. Time-distance 
techniques are used to study the characteristics of the transverse oscillations. We detect over 40 transverse 
oscillations in both bright and dark mottles, with periods ranging from 70 to 280 s, with the most frequent 
occurrence at ~ 165 s. The velocity amplitudes and transverse displacements exhibit characteristics similar 
to limb spicules. Neighbouring mottles oscillating in-phase are also observed. The transverse oscillations of 
individual mottles are interpreted in terms of magnetohydrodynamic kink waves. Their estimated periods and 
damping times are consistent with phase mixing and resonant mode conversion. 

Subject headings: Waves — magnetohydrodynamics (MHD) — Magnetic fields — Sun: Atmosphere — Sun: 
chromosphere — Sun: oscillations 



1. INTRODUCTION 

Dark and bright mottles are commonly observed in the 
quiescent solar chromosphere. They are small-scale jet- 
like features of relatively cool and dense material located 
at the boundaries of supergranular cells, ejected from the 
lower chromosphere at speed s of about 10-40 kms~' 
jTsiropo ula & Tziotzioul l2004t iRouppe van der Voort et alj 
|2007). A number of earlier studies suggest that bright and 
dark mottles correspond to phenomena in the lower and upper 
chromosphere, respectively ( Bray 1969), while others believe 
that the difference in brightness indicates different parts of the 
same structure (Beckers 1968, 1972; Banos & Macris 1970; 
[Bray & Loughheadlil974t ISterUng et al.lil993l) . Mottles are 
often con sidered as the disk rep r esentation of chromospheric 
spicu les (iHansteen et al.1 120061 : iRouppe van der Voort et alJ 
l2009h . The exact nature of mottles remains the sub- 
ject of an ongoing debate, with the majority of solar re- 
searchers agreeing that mottles and spicules are related, in 
the sense that they have simila r temperatures, dens ity profiles, 
widths, lengths and Ufe times (iTsiropoula & Schmi eder 1997]; 
IZachariadis et alJl200 1^. Chromospheric small-scale jet-like 
structures, such as mottles, spicules and fibrils, play an im- 
portant role in the mass balance of the solar atmosphere. It 
is estimated that only a small percentage of the mass outflow 
provided by mottles is sufficient to c ompensate for the coro- 
nal m ass loss due to the solar wind ( ITsiropoula & Tziotzioul 
l2004h . 

The highly dynamic photosphere can excite magnetohy- 
drodynamic (MHD) wave s which can p ropagate into the 
chromosphere and corona (Erdelvi 2006a,b). Spicular struc- 
tures can act as conduits for transferring wave energy from 
the lower to the upper parts of the solar atmosphere. Numer- 
ous observations of transverse motions in spicular structures 
have been reported in recent years (Kukhianidze et al. 200^ 



iDe Pontieu et'aD '2007'; 'He et al.' '2009a'; 'Tavabi et al.' "Mil 
Verth et aL 20 11; Okamoto & de Pontieu 201 1). An excellent 
review on the topic is given by Zaqarashvili & Erdelvi (2009). 
Kink waves have also been recently observed in active region 
fibrils (jPietarila et al.ll20Tll) . These are often considered as 
on disk, active region spicules which connect areas of op- 
posite magnetic polarity (van Noort & Rouppe van der Vo(3 
|2p06- 'De Pontieu et^ |2007|; |GuglielminoetaL| qOloT 
^uridze et al. 2011). Spicules and mottles are typically 
short-lived features, with lifetimes of the order of typical 
wave periods. The observed waves also appear to propagate 
with phase speeds close to the local Alfven speed. Mottles 
and spicules have typical lengths of 4 - 6 Mm, so waves 
propagate along the mottle in only 60 - 90 s. Both of these 
features make it challenging to obtain information on the ob- 
served wave motion, e.g., amplitudes, periods, phase speeds. 
However, the reported wave periods in spicules and fibrils are 
in the range 1-5 mins. Transverse motions in spicules have, 
in the past, been misinterpreted as Alfven waves, yet spicules 
(and mottles) appear as overdense structures in observations 
that are assumed to outline the quiet Sun magnetic field. 
In this case, the pure Alfve n waves are torsional motions 
(le.g.. Erdelvi & FedunI l2007h and the transverse waves are 
the MHD kink modes. 

Here we present high spatial and temporal resolution obser- 
vations of transverse oscillations in on-disk bright and dark 
mottles in the quiet solar chromosphere. We focus our at- 
tention on relatively long-lived mottles that display transverse 
motions which last for at least one wave cycle. We identify 
over 40 such waves and are able to obtain statistical informa- 
tion on the wave properties (e.g., amplitudes, periods). 



2. OBSERVATIONS AND DATA REDUCTION 
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Figure 1. Simultaneous ROSA images in Ha core (left), Ca ii K core (middle), and G band (right). Arrows indicate typical bright and dark mottles in Ha. White 
boxes indicate areas in which long-lived mottles have been detected in Ha. Axes are in heliocentric coordinates, where 1" « 725 km. 



The observations were obtained between 15:41 - 16:51 UT 
on 29 September 2010 with the Rapid Oscillations in the 
Solar Atmosphere (ROS A: Jess et a l. 2010) imaging system, 
mounted on the Dunn Solar Telescope (DST) at the National 
Solar Observatory, New Mexico, USA. A typical dataset in- 
cludes simultaneous imaging in G band. Ha core, and Ca ii K. 
High-order adaptive optics were used throughout the observa- 
tio ns. The ima ges were reconstructed using the algorithms 
of IWoger et alj {2008) followed by de-stretching. These al- 
gorithms were implemented to remove the effects of atmo- 
spheric distortion from the data. The effective cadence after 
reconstruction was 7.7 s for Ha, 9.6 s for Ca ii K, and 1 s for 
the G-band. The total field of view is 48" x 52" with a spatial 
sampling of 0.069"/pixel corresponding to a spatial resolution 
of 150 km in Ha. 

3. ANALYSIS AND DISCUSSION 

Figure [T| shows co-spatial and co-temporal images of the 
field-of-view in the Ha core, Ca II K and G-band. The Ho- im- 
age is dominated by elongated dark and bright mottles which 
appear co-spatial with Ca ii K brightenings and photospheric 
magnetic bright points. 

We find a total of 42 mottles (23 bright, 19 dark) that show 
transverse waves lasting for at least one wave cycle. In 7 
cases, neighbouring mottles are seen to oscillate in phase. 
Time-distance analysis reveals transverse oscillatory motions 
perpendicular to the mottle axis. The centre is determined 
by fitting a Gaussian profile to the cross-sectional flux pro- 
file for each time frame of the transverse cut (blue lines in 
the Figure |2]). The overall motion of the mottle shows a lin- 
ear trend which has been subtracted from the original sig- 
nal. Oscillation periods were determined using wavelet anal- 
ysis (Ijorrence & Compo 1998), and found to vary between 
70 - 280 s, with a strong peak at ~ 165 s. Transverse veloci- 
ties are in the range 3-18 km s"', with displacement ampli- 
tudes between 100 - 400 km (Figure|4]i. 

The phase speed of the transverse motions can be eval- 
uated by determining the phase difference between signals 
of the waves obtained at different heights along the mottle. 
Phase difference analysis along the mottle axis can be un- 
dertaken accurately only for 7 mottles due to the complex 
fine structure of the Ha images. Estimated phase speeds 
are in the range 40 - 110 km s ' . In most cases the oscilla- 
tions are seen only in a segment of the mottle typically less 
than 1 Mm long. A cadence of ~ 7.7 s limits the maxi- 



mum phase speed that may be detected to ~130 kms"'. The 
mean period, amplitude, and phase speed of the observed 
oscillations are similar t o results found for Umb spicules 
dKukhianidze et al.ll2006^ 'ZaqarashviU et al' 1 20071: iKim et all 
2008; He et al. 2009a; Okamoto & de Pontieii!|201lF 

Mottles can be considered as straight, cylindrical, high den- 
sity magnetic flux tubes which can guide different types of 
MHD wave modes, e. g., slow, fast, and Alfven. Alfven 
waves, propagating incompressibly along a magnetic flux 
tube, are torsional and do not displace the tube axis. MHD 
sausage modes may lead to intensity variations or variations 
of the mottle cross-section but do not to cause a transverse 
displacement of the tube axis. The clear transverse displace- 
ment of the mottle axis seen in our observations, suggests that 
the most appropriate interpretation for the oscillations is that 
of propagating/standing kink waves. The characteristic speed 
for the transverse mot ions is the kink speed, q, defined as 
(lEdwin & Robertsl[T983l) : 



Pi) + Pe 



(1) 



where po, vaq, Pe and v^e are the plasma densities and Alfven 
speeds inside and outside the tube, respectively. This equa- 
tion is applicable near the thin tube limit, i. e. r « A, where A 
and r are the wavelength and mottle radius, respectively. The 
wavelength. A, of the transverse wave is given hy A - Tc^, 
where T is the period of the wave. Hence A - 4800 km for 
r = 80 s and q - 60 km s"'. The radius of a mottle is typ- 
ically ~ 200 km, hence r A, the thin tube limit applies, 
and phase speed can be approximated by the kink speed. Fur- 
thermore, mottles are cool and dense material in the vicinity 
of strong magnetic flux concentration in network regions co- 
spatial with photospheric magnetic bright points (Figure [T]i. 
They may therefore, be considered as low-/? plasma magnetic 
tubes embedded in the field free (non magnetic) or weak field 
(compared to the mottle itself) environment. For no magnetic 
field outside the tube equation ([TJ, implies that 
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where Bq is the tube magnetic field strength. For a typ- 
ical value of the tube magnetic field Bo ~ 10 Gauss 
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Figure 2. Time-distance diagrams of the mottle intensity in Ho:. The blue line is the centroid of a Gaussian fit to the cross-sectional flux profile of the mottle at 
each time-step smoothed by a ~ 50 km width (red line). 
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Figure 3. The observed mottle displacements (black line) fitted with an exponentially damped cosine function (red line). 



jTruiill o Bueno et al J 120051). density po ^ (1 - 3) x 10"'^ 
g cm~ ^ (lBeckersll972l:ISterlingl2000l: lTsiropoula & Tziotziou^ 
120041) . and Pe/po ~ 0.1 - 0.033 we estimate from equation ^ 
kink speed range of ~ 46 - 86 km s This value compares 
well with our measured phase speeds (~ 40 - llOkms"'). 
The observed periods show a strong peak at around 165 s 
(top right panel of the Figure HJ. This value is very close 
to the chromospheric kink wave cut-off period defined as 
Tk = SHtt/q, where H is the chromospheric density scale 
heigh t (R ae& R oberts 1982; Spruit & Roberts 1983; Robert^ 
120041) . For the kink speed discussed above and a chromo- 
spheric scale height of H ~ 500 km, we estimate a kink 
wave cut-off period of ~ 170 s. Waves with periods 
above the cut-off become evanescent. However, wave prop- 
agation in non-ideal (e. g. with ion-neutral collisions or 
non-adiabatic effects) or inclined magnetic tubes can lead to 
an increase in the cut-off period and allows the penetration 
of low er-frequency oscillation s higher up in the solar atmo- 
sphere jDe Pontieu et al.f2004l) . Our imaging observations do 
not allow us to define accurate inclination angles for individ- 
ual mottles, however, as we are looking at the projection of 



mottle lengths onto the plane of the sky, this impUes that they 
are inclined from the vertical. Detection of periods above the 
cut-off value (Figure HJ suggests that these waves may tunnel 
through the evanescent region by propagating across inclined 
magnetic field lines. 

The damping properties of the observed transverse os- 
cillations are investigated by fitting the intensity profiles 
with an exponentially decaying cosine function, z(t) = 
Aexp(-f/T)cos(27rf/r + 0), where the fitting parameters 
A, T, T, (f) are the displacement am plitude, oscillation perio d 
and damping time, respectively (lAschwanden et all 120021) . 
We have found 12 cases where the mottle oscillations show 
evidence for damping. In Figure |3] we highlight 6 examples 
fitted with the above function and the oscillation parameters 
determined are listed in Table[T] The time sequence is divided 
into 5 shorter overlapping sections, and errors are obtained 
from the maximum and the minimum decay rates and phase 
of the sections. 

In Figure |5] we plot the estimated damping time as a func- 
tion of the oscillation period. We also calculate a best-fit 
power law scaling. The fit indicates an almost linear depen- 
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Figure 4. The lifetime (top left), period (top right), transverse velocity amplitude (bottom left), and transverse displacement amplitude (bottom right) for the 42 
mottles analysed in this study. Median values with standard deviations are also given. 



dence on the damping time of the kink wave period (Fig- 
ure |5]l. Resonant mode conversion predicts that the damping 
time is a hnear function of the period and may therefore be re- 
sponsible for the damping of the oscillations dGoossens et alj 
I200I iTerradas et aLll20ia ISoler et al.ll20Tlh in the chromo- 
sphere. Phase mixing follows a scaling law of t ~ T^^* (Fig- 
ure Is) that can also be a plausible the ory for the interpreta- 
tion of the observed damping (.Ofman & Aschwandenll2002t 
iMendoza-Briceflo et alj|2004l) . In contrast, a number of the 
transverse waves show an increase in ampUtude with time. 
This suggests the possibility of a co ntinuous driver or the in- 
fluence of a time -dependent plasma jMorton & Erdelyill2009l : 
lRudermanll201 II) . At present it is not possible to determine 
which, if any, of these options is responsible for the increas- 
ing amplitude. 

4. CONCLUDING REMARKS 

We have presented observations of relatively long-lived 
dark and bright mottles displaying transverse motions lasting 
for at least one wave cycle. The mottle oscillations studied 
in this work are interpreted as fast kink MHD waves. We ob- 
tain statistics on the properties, i.e. period, displacement and 
velocity amplitudes, of 42 separate events. Coherent oscil- 
lations in neighbouring mottles indicate an oscillatory source 
that operates on spatial scales larger than the mottle itself. The 



Table 1 

The amplitude A, oscillation period T, damping time t and phase for the 
damped oscillations. 



Damping 


A 


T 


T 




cases 


(km) 


(sec) 


(sec) 


(rads) 


1 


400±75 


165±12 


285±70 


0.45±0.15 


2 


150±50 


80±5 


62± 40 


0.42±0.22 


3 


300± 65 


190±5 


294±30 


0.54±0.1 


4 


330± 80 


170±17 


263±75 


2.8± 0.05 


5 


170± 60 


120±12 


191 ±35 


3±0.5 


6 


150± 50 


80±6 


266± 100 


3±0.1 


7 


170± 40 


80±20 


83± 25 


3±0.1 


8 


250 ± 50 


165 ±16 


476± 100 


2.9± 0.1 


9 


300± 60 


165±27 


136± 30 


3.4± 0.15 


10 


250± 70 


165±26 


200± 65 


0.12± 0.05 


11 


150±60 


170±14 


147± 65 


0.11 ±0.04 


12 


150±50 


260±10 


775 ± 170 


3.7± 0.3 



distance between in-phase oscillating mottles varies between 
300 - 1000 km, which is very similar to granular sizes. The 
relation between the damping times and periods suggests that 
resonant mode conversion and phase mixing may be viable 
damping mechanisms (Figure|5]l. Dissipation of the wave en- 
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Figure 5. Damping time r vs period T (red diamonds). Tlie .solid line is the 
best fit scaling law. The dashed line is the fit to r = ciT (scaling for resonant 
absoiption) with c\ = 2.27 and dotted line is the fit to t = C27'^^'* (scaling for 

erey asso transverse motions is thought to be 

important for the heating of the quiet solar coron a and the ac- 
celeration of the solar wind (le.g.. Cranmeiil200"7l) . 

Possible excitation mechanisms for the waves include 
granular buffeting, global oscillations and reconnec- 
tion events (Roberts 1979; Scruit 1981; Hollwee' '1981; 
Hasan & Kalkofen 1999 ; De Ponti eu et al., 2004; He et al.. 
2009bl) . Recently. iJesset al. (2012) suggested that longitudi- 



nal pressure modes in photospheric magnetic bright points 
can be converted into the transverse oscillations observed in 
Type I spicules. Fol lowing the suc cessful approa ch of coro- 
nal seismology (for reviews, see, e.g.. Ruderma n & Erdelyi l 



l2009t IVan Doorsselaere et alj 20091) . it may be possible that 



the observed kink waves could be exploited for seismology 
of the chromosphere (e.g., Verth et al. 2011). The higher 
temporal and spatial resolution of ground-based observing 
instruments (e.g., ROSA/DST, Crisp Imaging SpectroPo- 
larimeter (CRJSP)/Swedish Solar Telescope) provides a 
unique opportunity to study these short-lived, propagating 
chromospheric waves in detail. Such a study is not possible , 
at present, for the coronal counterpart dTomczyk et al.ll2007l ; 
lErdelvi & Ta royan 2008). 
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NSF, Hungary (OTKA, Ref. No. K83133). This work is sup- 
ported by the UK Science and Technology Facilities Council 
(STFC), with DBJ particularly grateful for the award of an 
STFC post-doctoral fellowship. We thank the Air Force Office 
of Scientific Research, Air Force Material Command, USAF 
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